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Abstract 

In  this  paper,  the  influence  of  gusts  on  the  dynamics  of  a 
large  flexible  aircraft  is  analyzed,  and  an  integrated  flight 
and  aeroelastic  control  law  that  reduces  gust  sensitivity  is 
presented.  The  calculations  are  based  on  an  integrated 
model  that  includes  all  1st  order  couplings  between  flight 
mechanic  and  structural  degrees  of  freedom.  Uniform,  1- 
dimensional  and  multidimensional  gust  models  are  imple- 
mented and  used  for  gust  sensitivity  analysis.  For  the  exam- 
ple aircraft,  the  differences  in  gust  sensitivity  calculated 
with  the  1-dimensional  and  multi-dimensional  gust  models 
are  significant  Integrated  attitude,  stability  augmentation, 
and  aeroelastic  control  laws  for  longitudinal  and  lateral 
motion  are  designed  using  (i  - synthesis.  With  the  control 
laws,  flight  maneuvers  do  not  excite  elastic  reactions,  and 
the  sensitivity  to  gusts  is  considerably  reduced. 

Nomenclature 

Bhh  = modal  damping  matrix 

k = reduced  frequency 

Khh  = modal  stiffness  matrix 

L = characteristic  wave  length  for  gust  (2500  ft) 

Mh) , = modal  mass  matrix 

Ph  = modal  applied  aerodynamic  forces  (gust) 

Q(ik)  = modal  aerodynamic  force  coefficient  matrix 
tabulated  for  reduced  frequencies  k 
q = state  vector  in  generalized  coordinates 

xj  = x location  of  aerodynamic  panel 

x0  = reference  value  for  aerodynamic  coordinate  system 
«,v,w  = longitudinal,  lateral  and  vertical  velocity 

p,q,r  = roll,  pitch  and  yaw  rate 

w = aerodynamic  downwash 

a , |3 , y = angle  of  attack,  sideslip  angle  and  flight  path  angle 

Y = dihedral  angle 

<t>  = mode  shape  matrix 

<j> , 0 , V|/  = roll  attitude,  pitch  attitude  and  heading 

<p  = phase  angle 

d>  = matrix  of  one  dimensional  spectrum  function 

T*  = matrix  of  two  dimensional  spectrum  function 

(U  ,^so  = symmetric  and  antimetric  deflection  of 

^ai  ’ ^ao  inner  and  outer  ailerons 

i;  ,ri , C,  = aileron,  elevator  and  rudder  deflection 
to  = frequency 

Q.  = co/V  wave  number 

1 Introduction 

Today’s  airlines  are  requesting  bigger  and  more  fuel  efficient 

aircraft  to  reduce  their  operation  costs.  Consequently,  fuselages 

and  spans  of  new  aircraft  designs  are  getting  longer,  and  the 

need  to  reduce  structural  weight  reduces  structural  stiffness. 


Both  effects  lead  to  more  flexible  aircraft  structures  with  sig- 
nificant aeroelastic  coupling  between  flight  mechanics  and 
structural  dynamics,  especially  at  high  speed,  high  altitude 
cruise.  To  counteract  gust  and  maneuver  induced  aeroelastic 
vibrations,  which  impair  ride  comfort  and  structural  loads, 
active  aeroelastic  control  strategies  are  investigated. 

Using  separate  models  for  flight  mechanics  and  aeroelastics 
(sufficient  for  the  smaller  and  more  rigid  aircraft  in  service 
today)  aeroelastic  coupling  cannot  be  described.  Further,  it  is 
doubtful  whether  simple  uniform  or  ID  gust  models  are  still 
adequate  for  large  aircraft.  To  address  the  former  deficiency, 
integrated  models  describing  both  the  flight  mechanic  and  the 
aeroelastic  behavior  of  flexible  aircraft  have  been  developed 
recently1.  The  issues  of  gust  modelling  and  active  aeroelastic 
control  are  addressed  here. 

Calculations  of  gust  loads  are  of  major  importance  for  the 
structural  and  aerodynamic  layout  of  an  airplane.  Certification 
requirements  due  to  discrete  gusts  as  well  as  continuous  gust 
loads  have  been  developed.  FAR  252  describes  one  dimen- 
sional (ID)  1-cos  shape  gust  requirements.  The  FAA-ADS  532 
report  describes  ID  continuous  gust  requirements,  which  are 
based  on  power  spectral  density  PSD  methods. 

There  always  have  been  efforts  to  extend  existing  gust  formu- 
lations to  a more  realistic  modelling  as  the  most  widely  used 
ID  gusts  are  the  most  abstract  realization  of  real  gusts,  which 
will  be  random  not  only  in  one  direction  but  in  each  direction. 
The  energy  distribution  within  gusts  has  already  been 
described  in  the  early  60s  by  Dryden  and  von  Karman.  These 
gust  spectra  have  been  developed  in  3D  space,  but  in  the  later 
years  only  the  ID  spectra  have  been  used.  Nevertheless, 
research  on  multidimensional  formulations  of  gusts  and  the 
modelling  of  the  induced  aerodynamic  forces  is  continuing.  On 
the  evaluation  of  multidimensional  gust  fields  and  the  resulting 
airplane  responses,  most  authors  rely  on  cross-spectral  formu- 
lations. Crimaldi3  describes  a method  for  2D  gust  modelling 
which  is  based  on  calculations  of  the  cross-spectrum  of  a gust 
acting  on  single  aerodynamic  strips  of  panels  along  the  longitu- 
dinal axis  of  the  airplane.  These  gust  cross-spectra  depending 
on  the  lateral  separation  distance  have  been  developed  by 
Eichenbaum4  and  are  described  by  Bessel  functions  of  the  third 
kind. 

The  gust  models  described  in  this  paper  avoid  the  evaluation  of 
cross-spectral  density  functions  by  applying  directly  ID  and 
2D  PSD  functions  to  vertical  and  lateral  gusts  as  in  Etkin5. 
Starting  from  a very  basic  formulation  gust  models  are  step  by 
step  developed  towards  a fully  multidimensional  gust  descrip- 
tion, so  that  effects  caused  by  the  increase  in  accuracy  of  the 
models  can  be  identified. 

Conventional  flight  control  systems  are  designed  to  assure 
good  handling  qualities  and  to  eliminate  the  influence  of  elas- 
ticity, treating  the  aircraft  as  a rigid  body.  For  ride  comfort 
improvement,  separate  aeroelastic  control  loops  have  been 
implemented  on  some  modem  aircraft6.  The  flight  control  and 
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the  aeroelastic  control  loops  are  then  separated  by  dynamic  fil- 
ters. As  rigid  body  dynamics  and  low  frequency  elastic  modes 
get  closer  with  increasing  structural  flexibility,  separate  design 
of  stability  augmentation  systems  and  aeroelastic  control  loops 
becomes  more  difficult.  Therefore,  several  recent 
studies1,7’8,9,10  have  investigated  integrated  flight  mechanic  and 
aeroelastic  control  design.  In  this  study,  an  integrated  flight  and 
aeroelastic  control  law  for  a heavy  study  four-engine  civil 
transport  aircraft  at  high-speed,  high-altitude  cruise  is  pre- 
sented, and  the  influence  of  the  integrated  control  law  on  gust 
sensitivity  is  analyzed. 

2 Integrated  aircraft  model 

Linearized  integrated  flight  mechanic  and  aeroelastic  models 
for  different  flight  (20000ft  to  30000ft  altitude  and  Ma=0.4  to 
Ma=0.86)  and  load  (full  and  empty  wing  and  trim  tank)  condi- 
tions are  derived  for  the  example  aircraft.  The  modelling  fol- 
lows the  procedure  described  by  Schuler1. 

Structural  Dynamics 

The  structural  dynamics  of  the  aircraft  structure  are  analyzed 
using  a detailed  full-span  FEM-model.  A normal  mode  analy- 
sis is  performed  and  the  low-frequency  elastic  modes  up  to  a 
frequency  of  about  20Hz  are  retained.  Additional  mode  shapes 
for  rigid  body  motion  and  control  surface  deflections  are  gen- 
erated. Modal  coordinates  q(qr  for  rigid  body  motion  and  qe 
for  elastic  modes)  are  introduced  and  used  to  describe  rigid 
body  and  structural  dynamics: 


domain.  For  this  study,  the  Minimum  State  Method  of  Karpel14 
has  been  used  to  perform  the  approximation. 


Fig.  1 : Aerodynamic  model  (NASTRAN)  - the  lifting  surface 
is  discretized  using  1570  doublet  elements  (boxes) 

For  the  time-domain  model,  gusts  are  modelled  as  acting  uni- 
formly on  the  aircraft  structure  and  separately  on  longitudinal 
and  lateral  motion.  Then  a (linear,  time-invariant)  state-space 
description  of  the  (linearized)  aircraft  dynamics 


x = Ax  + Bu  + Ew 
z = Cx  + Du  + Fw 


(4) 


with  control  inputs  u and  (uniform)  gust  inputs  w is 
derived12,13,1.  The  uniform  gust  model  however  is  too  inaccu- 
rate an  approximation  for  gust  analysis.  Therefore,  in  the 
sequence,  more  realistic  gust  models  are  discussed  (sections  4 
and  5).  Finally  it  is  analyzed  how  the  controllers  designed  with 
the  uniform  gust  model  perform  with  respect  to  more  realistic 
gust  models  (section  9). 


ir=f(ir)  + E?  + F?  + FA  (1) 

he  = -2CmA<?e  - A 2qe  + FF  + FA 

with  A = diag(co,) , modal  damping  coefficient  , (here 
= 0.01 ),  aerodynamic  forces  FA , thrust  forces  FF  and 
gravitational  forces  Fa . Displacement,  velocity  and  accelera- 
tion of  the  structure  in  physical  degrees  of  freedom  (z,  z,  z ) are 
expressed  in  terms  of  modal  coordinates  by 

Z = <Pq  , z = <i>q , z = dxq , (2) 

where  <J>  is  the  mode  shape  matrix. 

Aerodynamics 

Steady  and  unsteady  aerodynamic  forces  are  calculated  using 
the  Doublet-Lattice  method11,  as  implemented  in  NASTRAN. 
The  location  of  the  1 570  boxes  (each  representing  a doublet- 
singularity)  is  shown  in  Fig.  1 . Aerodynamic  force  coefficients 
are  calculated  for  18  discrete  reduced  frequencies  and  tabu- 
lated in  matrices  Q{ik) . Linearized  aerodynamic  forces  can 
then  be  written  as 


Fr 

F* 


Qrr(ik)  Qre{ik) 

Qr 

+ 

Qrgm 

Qer(ik)  Qee{ik)_ 

.q1 

Qeg(ik) 

with  Qre  and  Qer  representing  the  cross-coupling  between 
rigid-body  and  structural  modes  and  Qrg  and  Qeg  the  forces 
due  to  gust  vector  inputs  w . In  the  form  of  Eq.  3,  aerody- 
namic forces  can  only  be  evaluated  for  harmonic  oscillations  at 
given  discrete  frequencies.  For  simulation  purposes  and  control 
law  design,  the  aerodynamic  force  description  is  therefore 
extended  to  the  Laplace  domain,  approximated  by  a rational 
transfer  function  matrix  and  then  transformed  to  the  time 


3 Random  turbulence  models 

For  most  evaluations  turbulence  is  regarded  as  a random  phe- 
nomenon. Its  energy  distribution  has  to  be  described  in  the  3D 
space. 


Fig.  2:  2D  spectrum  ¥33  of  vertical  von  Karman  turbulence 

Considering  vertical  gusts  for  airplane,  it  is  apparent  that  only 
variations  in  longitudinal  and  lateral  direction  will  be  signifi- 
cant and  variations  of  gust  velocity  in  vertical  direction  can  be 
neglected  as  the  airplane  height  is  small  in  comparison  to  span 
and  length.  2D  spectra  will  therefore  be  appropriate.  In  this 
study,  the  von  Kdrman  Spectrum15  is  used  to  describe  the 
energy  distribution  versus  frequency.  For  vertical  turbulence  it 
is  given  by 


¥ 33(^1*  ^2)  — 


16c2(aL)4 

9n 


(0\  + Q22) 


2 2 2 2 7/3 

(1  +aT(fl;  + Oj)) 


• (5) 


Results  presented  in  this  paper  will  be  calculated  with  the 
appropriate  ID  and  2D  von  Karmdn  spectra.  Fig.  2 shows  the 
2D  von  KSrmSn  spectrum.  ID  spectra  are  obtained  by  integrat- 
ing Eq.  5 with  respect  to  Cl2-  Wavenumber  Q2  ar>d  £2 ; will  be 
discussed  in  section  5. 
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4 Standard  ID  continuous  gusts  response 
For  ID  gust  loads  only  changes  of  gust  velocity  in  flight  direc- 
tion are  considered,  for  vertical  as  well  as  for  lateral  gusts. 
One  way  to  compute  these  gust  forces  is  to  assume  a uniform 
gust  velocity  over  the  airplane.  Gust  forces  are  obtained  that 
are  equal  to  the  generalized  aerodynamic  forces  of  the  transla- 
tional vertical  and  lateral  rigid  body  modes.  Standard  gust  cal- 
culations (NASTRAN)  consider  phase  delay  effects  along  the 
flight  path.  The  gust  field  is  assumed  to  be  frozen  (Taylor’s 
hypothesis).  Gust  forces  are  related  to  a normal  downwash  at  a 
specific  aerodynamic  box  j by 

Pk(to)  = Qhj(Ma,k)Wj((o),  (6) 

where  Qkj  relates  downwash  and  modal  aerodynamics  forces 
and  is  available  through  NASTRAN.  The  downwash  of  a verti- 
cal gust  which  actually  describes  the  phase  delay  is  given  by16 


w;((0)  = cos  , (7) 

where  the  normal  downwash  on  a aerodynamic  box  is  propor- 
tianal  to  the  dihedral  angle  y,-.  For  lateral  gusts  cosyj  has  to  be 
replaced  by  siny^. . 


Fig.  3:  PSD  of  acceleration  of  outer  wing  due  to  uniform  ver- 
tical gust  and  gust  with  phase  delay  effects  in  flight 
path  direction 

The  flutter  equation  in  generalized  coordinates  is  used  for  open 
loop  simulations  with  applied  forces: 

[-Mhha2  + iBhhi o + Khh  - gp  V^Q^Ma,  *)]?  = Ph{ coX8) 

In  this  equation  modal  masses  Mkk,  damping  Bhh  and  stiffness 
Bhk  are  real  matrices,  the  matrix  of  modal  aerodynamic  forces 
Qhk  coefficients  is  complex.  The  result  is  obtained  in  modal 
coordinates.  Displacement,  velocity  and  acceleration  of  the 
structure  in  physical  degrees  of  freedom  can  be  computed  with 
Eq.  2.. 


Fig.  4:  PSD  of  acceleration  of  fuselage  close  to  pilot  and  rear 
fuselage  due  to  ID  lateral  gust  with  phase  delay 

The  RMS  values  of  the  response  in  physical  coordinates  are 
obtained  by  applying  ID  gust  spectra.  Considering  the  phase 


delay  along  the  flight  path  is  the  first  step  towards  a 2D  gust 
model.  The  effect  on  the  acceleration  of  the  wing  can  be  seen 
in  Fig.  3.  The  evaluation  of  the  response  of  the  wing  bending  is 
used  in  most  studies  for  the  evaluation  of  the  critical  loads  of 
an  airplane.  This  might  be  in  general  true  for  smaller  aircraft, 
but  from  Fig.  3 two  important  things  can  be  seen:  Phase  delay 
changes  the  aircraft  response,  and  there  are  modes  beside  sym- 
metric wing  bending  at  higher  frequencies  that  will  be  consid- 
erably excited  independent  of  gust  model. 

As  to  the  response  of  an  airplane  to  lateral  gusts,  it  is  obvious 
that  the  so  called  “fish  tailing”  mode  and  “fuselage  bending” 
mode  at  higher  frequency  should  be  considered.  The  “fish  tail- 
ing” mode  results  in  high  accelerations  in  the  rear  fuselage  of 
the  airplane  and  small  loads  at  the  front  fuselage,  see  Fig.  4. 
For  clarity:  response  of  uniform  gust  is  not  depicted. 


5 Multidimensional  continuous  gust  response 
Gust  calculations  for  multidimensional  continuous  gusts  are 
based  on  the  method  described  by  Etkin5.  It  differs  from  the  ID 
continuous  gust  in  that  there  is  an  additional  variation  along  the 
span.  The  gust  velocity  is  varying  along  the  flight  direction  as 
before  and  is  “wavy”  over  the  airplane. 


Fig.  5:  Aircraft  flying  through  a 2D  vertical  gust  field 

For  the  vertical  gust  case  a 2D  gust  field  is  obtained,  where  the 
downwash  is  dependent  on  wave  numbers  Slj=a/V  and  Q2. 
Eq.  7 is  modified: 


(£5,(x,-j:0)  + £1  2y,} 

= cosYj-e 


(9) 


It  can  be  seen  that  the  second  wave  number  Q2  results  in  a dif- 
ferent phase  of  the  downwash  compared  to  the  ID  gust 
assumption.  Assuming  a fixed  proportion  of  £22  and  Q;  the  sit- 
uation can  also  be  described  by  the  airplane  flying  over  a ID 
gust  field  with  an  angle  a,  see  Fig.  5.  Lateral  gusts  are  imple- 
mented similar  to  vertical  gusts.  By  implementing  the  depen- 
dence on  n2  Eq.  8 is  extended  to  a 2D  transfer  function.  As  an 
example  the  acceleration  of  the  wing  for  different  values  of  £22 
is  presented.  It  can  be  shown  that  dependent  on  each  character- 
istic mode,  the  second  wave  number  Q2  yields  higher  or  lower 
accelerations.  Additionally,  it  has  to  be  mentioned  that  because 
of  the  phase  delay  the  response  of  the  left  and  right  side  of  the 
airplane  will  be  different.  Similar  to  the  lD-case  the  response 
has  been  multiplied  with  the  2D-spectrum.  Neglecting  longitu- 
dinal gusts  it  has  been  shown5  that  the  total  RMS  response  for 
isotropic  turbulence  is  given  by  the  uncorrelated  response  of 
vertical  and  lateral  gust.  Therefore  the  spectra  for  lateral  and 
vertical  motion  can  be  computed  separately  and  added  : 


V/Q>’n2)  = N Tw(ni'  n2> + |xiW!2'?KW(Qi-  n2> 


(10) 


In  this  equation  'Fxx  is  the  2D  spectrum  response  of  a selected 
physical  degree  of  freedom  xt.  Xiv  is  the  response  due  to  unit 
gust  input. 
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Frequency 


Fig.  6:  PSD  of  acceleration  of  outer  wing  due  to  2D  vertical 
gust  for  £22=0.02  and  Q2=0.08 


quency.  As  at  higher  frequencies  gust  spectra  become  very 
small,  only  the  first  number  of  elastic  modes  will  lead  to  PSD 
values  that  are  used  for  evaluation. 

Adding  lateral  gusts  to  the  response  according  to  Eq.  10  causes 
a slight  increase  of  loads  on  surfaces  that  are  mainly  affected 
by  vertical  gusts  (see  Fig.  8).  For  antimetric  modes  (“fish  tail- 
ing”) the  situation  will  be  vice  versa. 


Fig.  8:  PSD  of  vertical  acceleration  of  outer  wing  due  to  ver- 
tical 2D  gusts  alone  and  a combination  of  vertical  and 
lateral  2D  gusts 


To  compare  these  results  with  the  ID  case  we  need 

«(«,)  = J ^(n,,  n2)<«2.  (11) 

With  this  equation  the  acceleration  of  the  wing  and  of  the  outer 
engine  due  to  ID  and  2D  gust  input  is  compared  and  presented 
in  Fig.  7. 


Frequency 


Fig.  7:  PSD  of  vertical  wing  and  lateral  outer  engine  accelera- 
tion for  ID  and  2D  gust  input 

It  can  be  seen  that  the  2D  formulation  leads  to  lower  accelera- 
tions for  the  first  wing  bending  mode.  Due  to  the  antimetric 
phase  delay  antimetric  modes  are  excited,  which  can  clearly  be 
seen  for  the  first  antimetric  wing  bending  mode  . This  antimet- 
ric response  of  the  first  wing  bending  explains  the  reduction  in 
the  symmetric  wing  bending,  as  the  ID  spectrum  with  fi2=0 
puts  too  much  emphasis  on  the  first  wing  bending  mode.  In 
contrast  to  the  first  wing  bending  mode,  there  are  modes  whose 
response  can  increase  for  the  2D  gust  formulation.  The  reason 
is  that  based  on  the  phase  delay  there  are  situations  when  the 
gust  shape  over  the  wing  is  equal  to  a characteristic  wing  bend- 
ing form  and  will  therefore  amplify  the  response  at  this  fre- 


6 Open  loop  system  analysis 


With  all  models  based  on  NASTRAN  frequency  domain  aero- 
dynamics, frequency  responses  are  the  natural  choice  for  air- 
craft dynamics  analysis  and  offer  the  most  accurate  results. 
Further  advantages  are  the  possibility  to  clearly  identify  the 
influence  of  individual  elastic  modes,  to  assess  the  influence  of 
random  disturbances  (such  as  gusts)  and  to  readily  draw  con- 
clusions for  control  design  (stability  margins  etc.).  While  fre- 
quency responses  could  also  be  obtained  from  the  linear  state- 
space  models  (Eq.  4),  those  shown  below  have  been  calculated 
by  directly  solving  the  flutter  equation  (see  Eq.  8),  thus  avoid- 
ing errors  resulting  from  the  transformation  to  the  time- 
domain. 


Elevator  on  cockpit  vertical  acceleration  r|  -4  ipi 


Fig.  9:  Cockpit  vertical  acceleration  frequency  response  to 
elevator  input  for  different  flight  conditions 

Fig.  9 shows  the  transfer  function  from  elevator  to  cockpit  ver- 
tical acceleration  for  different  flight  conditions.  As  expected 
from  rigid-body  flight  mechanics,  the  gain  amplitude  increases 
with  dynamic  pressure,  the  phugoid  frequency  decreases,  and 
the  frequency  of  the  short  period  motion  increases  with  air 
speed.  In  the  aeroelastic  frequency  range,  the  gain  response 
shows  the  influence  of  dynamic  pressure  on  amplitude,  while 
phase  response  is  not  much  affected  by  changes  in  flight  condi- 
tion. For  changes  in  load  condition,  however,  the  situation  is 
different.  The  load  distribution  (fuel  and  payload)  strongly 
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influences  the  dynamic  behaviour  of  the  elastic  structure  and 
consequently  aeroelastic  coupling.  Fig.  10  shows  the  frequency 
response  for  three  different  load  conditions  at  constant  flight 
condition.  While  changes  in  elastic  mode  shapes  and  frequen- 
cies are  not  unexpected  (wing  bending  frequency  should 
increase  with  fuel  consumption),  strong  changes  in  damping 
(maximum  amplitudes)  and  phase  response  are  also  observed. 


Fig.  10:  Cockpit  vertical  acceleration  frequency  response  to 
elevator  input  for  different  load  conditions 

A fixed  gain  control  design  therefore  faces  a robustness  chal- 
lenge in  both  the  low-frequency  range  (rigid  body  dynamics 
changing  mainly  with  flight  condition)  and  the  high-frequency 
range  (aeroelastics  changing  mainly  with  load  condition).  With 
the  low-frequency  elastic  modes  being  close  to  the  short  period 
mode,  (the  1st  wing  bending  mode  is  not  visible  in  Figs.  9 and 
10)  strong  coupling  between  rigid  body  and  structural  dynam- 
ics is  observed. 

Under  these  premises,  the  aforementioned  difficulties  associ- 
ated with  separate  flight  and  aeroelastic  control  and  the  advan- 
tages offered  by  integrated  control  are  evident. 

7 Control  design  for  integrated  flight  and 
aeroelastic  control 

Actuators 

The  integrated  control  law  commands  the  conventional  control 
surfaces  for  primary  flight  control,  i.e.  elevator,  rudder,  and 
inner  and  outer  ailerons.  While  Schuler1  assumes  symmetri- 
cally deflected  inner  ailerons  to  be  available  in  the  longitudinal 
motion,  in  this  study  symmetric  inner  and  outer  aileron  deflec- 
tion is  made  available,  however  restricted  to  low  authority 
aeroelastic  control  purposes. 

Sensors 

As  all  physically  realizable  sensor  signals  contain  both  rigid 
body  and  elastic  motion,  sensor  number  and  position  is  an 
important  consideration.  Kubica  and  Livet8  and  Ward  and  Ly7 
use  the  signals  of  the  conventional  aircraft  sensor  platform.  In 
this  study,  roll  attitude  <|>  and  pitch  attitude  0 , roll,  pitch  and 
yaw  rate  (p,q,r ),  vertical  and  lateral  acceleration  are  assumed  to 
be  measured  by  a sensor  platform  in  the  forward  fuselage  sec- 
tion. As  in  Schuler1,  additional  accelerometers  at  different 
points  of  the  aircraft  structure  are  assumed  to  be  available. 
They  are  placed  at  positions  where  low-frequency  signals  are 
dominant  (providing  for  a physical  low-pass  filter)  using  a sen- 
sor placement  strategy17.  Vertical  acceleration  at  a mid-wing 


position  and  lateral  acceleration  at  inner  and  outer  engines  are 
retained  for  both  longitudinal  and  lateral  control.  Lateral  accel- 
eration at  a rear  fuselage  position  is  added  for  lateral  control. 

Model  reduction 

As  the  linearized  integrated  model  (see  section  2)  features  a 
high  number  of  states  (e.g.  103  states  for  the  longitudinal 
motion)  an  order  reduction  is  required  to  derive  a control 
design  model.  A 30-state  reduced  order  model  is  obtained  by 
applying  a combination  of  balanced  truncation  and  balanced 
condensation  model  reduction  techniques18. 

Choice  of  the  Design  Method 

p-synthesis 19,20  has  been  chosen  for  control  design  for  several 
reasons.  Since  the  aircraft  model  is  based  on  frequency-domain 
aerodynamics,  frequency  domain  control  design  is  advanta- 
geous, as  important  time-domain  model  characteristics,  such  as 
pole  positions  do  not  directly  relate  to  the  physical  modelling 
process  but  result  from  the  transformation  of  the  frequency 
domain  model  to  the  time  domain.  The  interpretation  of  short 
period  frequency  and  damping  for  example  is  not  possible  in 
terms  of  short-period-pole-position  alone,  because  the  influ- 
ence of  aerodynamic  lag-states  has  to  be  considered. 

A MIMO  control  design  method  has  been  selected,  as  com- 
bined control  surface  actuation  is  much  more  effective  for 
aeroelastic  control  than  separate  control  loops  for  each  control 
input,  given  the  position  of  elevator  and  ailerons  on  the  aircraft 
structure.  With  30-state  control  design  models,  the  availability 
of  numerically  reliable  algorithms  for  Hx  - optimization  / p- 
synthesis  is  another  important  consideration.  Closed  loop 
shaping21,22  has  been  chosen  for  setting  up  the  HM  - optimiza- 
tion problem  for  its  flexibility  in  translating  design  specifica- 
tions into  weighting  functions. 

The  design  process  follows  the  procedure  outlined  in  Ref.  10. 
The  C*  combination  of  pitch  rate  and  vertical  acceleration 
(C*  = q-z/V)  is  chosen  as  control  variable  for  the  longitu- 
dinal motion.  A proportional-integral  control  characteristic  is 
achieved  by  adding  the  integral  of  C*  , C*  as  an  additional 
measurement.  Accordingly,  in  the  lateral  motion,  roll  rate  p is 
the  control  variable  and  roll  angle  <(>  is  added  for  Pi-control. 

Command  tracking,  vibration  reduction,  aeroelastic  damping 
increase,  control  effort,  disturbance  attenuation,  and  gust  sensi- 
tivity are  considered  for  nominal  performance  in  the  design  by 
including  and  weighting  the  corresponding  transfer  functions 
in  the  performance  index  for  HM-  control.  As  the  more  sophis- 
ticated gust  models  described  in  this  article  are  not  available 
for  (state-space-model-based)  control  design,  the  gust  sensitiv- 
ity objectives  are  formulated  using  the  uniform  gust  approach. 
In  a robust  performance  formulation,  robustness  is  demanded 
against  changes  in  the  aeroelastic  frequency  and  damping 
parameters  and  strengthened  in  several  p-synthesis  design  iter- 
ations. The  p -synthesis  design  is  executed  using  the  Xmath- 
Software  package. 

The  resulting  40th-order  longitudinal  controller  features  rela- 
tively high  low-frequency  gains  for  pitch  rate  and  C*  feedback 
to  elevator,  nearly  no  low-frequency  aileron  activity  and  a rap- 
idly decreasing  control  effort  in  the  frequency  range  of  the 
elastic  modes.  Acceleration  measurements  are  mainly  used  for 
aeroelastic  control  purposes  in  the  lower  frequency  range. 
Using  Hankel  norm  optimal  model  reduction23,  a controller 
order  reduction  to  26  states  can  be  performed  without  signifi- 
cant performance  losses.  High-frequency  control  gain  attenua- 
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tion  however  is  lowered.  If  the  controller  order  is  further 
reduced,  performance  rapidly  deteriorates.  In  the  lateral 
motion,  roll  rate  and  roll  angle  feedback  on  ailerons  and  yaw 
rate  and  sideslip  (estimated  from  lateral  acceleration)  feedback 
on  rudder  dominate  the  integrated  control  law  as  expected. 
Again,  the  feedback  of  the  acceleration  measurements  is  prin- 
cipally used  for  aeroelastic  control  in  the  low  aeroelastic  fre- 
quency range.  Model  reduction  techniques  are  applied  to 
reduce  the  controller  order  (to  26  states). 

8 Closed  loop  system  analysis 

In  this  section,  an  analysis  of  the  closed  loop  aircraft  dynamics 
is  presented.  The  discussion  is  intended  to  give  an  overview 
over  the  different  performance  objectives  considered  for  p- 
synthesis.  For  the  sake  of  brevity,  it  is  restricted  to  the  longitu- 
dinal motion.  The  results  achieved  with  the  lateral  control  law 
are  similar  to  previous  results10.  Robustness  to  system  changes, 
however,  has  been  increased  considerably  by  using  p-synthesis 
instead  of  standard  control. 

The  robust  integrated  flight  and  aeroelastic  controller  consider- 
ably increases  the  damping  of  the  low-frequency  elastic  modes 
as  shown  by  the  comparison  of  open  and  closed  loop  pole  loca- 
tions in  Fig.  11.  Damping  increase  is  40  to  400%  for  the  first  6 
elastic  modes. 
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Fig.  11:  Pole  shifting  - longitudinal  control;  final  p - controller 
on  linearized  plant  model,  Mach=0.86,  b= 30000ft; 
plant  pole  damping  is  increased,  while  compensator 
(observer)  poles  are  shifted  to  the  right 

For  geometric  reasons  it  could  be  expected  that  the  damping  of 
wing  and  engine  modes  is  mainly  due  to  the  availability  of 
symmetric  aileron  control  input,  whereas  the  damping  of  fuse- 
lage bending  modes  is  achieved  by  elevator  activity.  While  this 
was  generally  confirmed  in  the  design  process,  coordinated 
control  surface  deflection  enhances  performance.  For  flight 
and  load  conditions  beyond  the  design  case,  the  damping 
increase  achieved  with  the  robust  fixed  gain  controller  in  the 
loop  and  control  performance  in  general  decreases.  Although, 
in  the  present  case,  a damping  increase  can  be  achieved  for  all 
flight  and  load  conditions,  the  limits  of  robust  fixed  gain  con- 
trol for  changing  flight  conditions  are  evident. 


The  closed  loop  time  response  to  a Is  - pulse  command 
y=0.027s  and  a 3.5°  symmetric  aileron  disturbance  inpulse  at 
r=6s  is  shown  in  Fig.  12  (empty  symbols). . 
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Fig.  12:  Closed  loop  response  comparison  of  a flight  controller 
and  the  integrated  flight  and  aeroelastic  controller 

For  an  assessment  of  the  robust  integrated  controller,  the  time 
response  of  a simple  C*  -flight  controller  is  shown  in  compari- 
son (filled  symbols).  Command  tracking  is  equally  fast  for  both 
controllers,  with  a somewhat  lower  short  period  frequency  for 
the  integrated  controller. 
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Fig.  13:  Open  and  closed  loop  sensitivity  to  input  perturbations 
- p-controller  longitudinal  motion;  Singular  values  of 
G and  SyG;  Sensitivity  to  input  disturbances  is 
reduced  in  the  frequency  range  of  controller  operation; 

The  control  effort  of  the  integrated  control  law  is  based  mainly 
on  elevator.  Low  amplitude  inner  and  outer  aileron  deflections 
are  used  for  aeroelastic  damping  (but  not  for  direct  lift  control). 

While  command  tracking  excites  practically  no  elastic  vibra- 
tions, the  disturbance  impulse  causes  aeroelastic  vibrations.  It 
can  be  seen  that  the  integrated  controller  rapidly  damps  theses 
vibrations,  especially  in  the  fuselage,  but  also  at  the  inner 
engine.  The  handling  qualities  achieved  with  this  controller  are 
significantly  improved  over  previous  results10. 
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Input  sensitivity 

The  singular  values  of  the  transfer  matrix  from  the  control 
(internal)  inputs  to  the  measurement  (internal)  outputs  are  cal- 
culated for  both  the  open  loop  (G)  and  the  closed  loop  system 
( SyG ) to  analyze  the  sensitivity  to  input  perturbations. 

Fig.  13  shows  the  singular  value  diagram  for  the  design  flight 
and  load  condition  for  rigid  body  and  elastic  mode  frequencies. 
It  can  be  seen  that  the  input  disturbance  sensitivity  is  consider- 
ably reduced  in  the  frequency  range  of  the  phugoid  and  short- 
period  motion,  the  lightly  damped  engine  modes,  and  the  fuse- 
lage bending  mode.  The  input  sensitivity  in  the  range  of  higher 
frequency  elastic  modes  is  not  affected  by  the  control  law.  This 
observation  confirms  the  conclusions  drawn  from  the  pole- 
shifting  analysis. 

9 Closed  loop  gust  sensitivity 

In  this  section,  a detailed  analysis  of  closed  loop  performance 
with  respect  to  gusts  is  presented.  In  the  first  run,  the  uniform 
gust  models  used  for  control  design  and  the  NASTRAN  ID 
gust  model  (standard  for  smaller  aircraft)  are  employed.  In  a 
second  analysis  run,  the  2D  gust  model  of  section  5 is 
employed. 

Gust  sensitivity  analysis  using  ID  gust  models 

For  the  gust  sensitivity  analyses  of  the  first  run,  a ID  von 
Karman  gust  spectrum15  is  used. 


Fig.  14:  ID  vertical  gust  response  - vertical  acceleration  in  the 
cockpit;  open  and  closed  loop  sensitivity,  Mach  0.7, 
/i=25000ft,  p-controller;  the  gust  sensitivity  in  the  fre- 
quency range  of  the  rigid  body  motion  is  reduced 
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Fig.  15:  Open  and  closed  loop  sensitivity  to  vertical  gusts; 
a)  uniform  gust  model;  b)  ID  gust  model; 
von  Kdrmdn  spectrum;  p-controller  longitudinal 
motion,  Ma=0.7,  A=25000ft;  the  gust  sensitivity 
reduction  in  the  range  of  engine  and  fuselage  modes  is 
significant,  independently  of  the  gust  model  employed 

The  same  analysis  is  performed  for  the  lateral  control  law.  Fig. 
16  shows  the  cockpit  lateral  acceleration  response  to  a von 
KSrmdn  spectrum  - based  ID  lateral  gust  input  for  the  open  and 
the  closed-loop  system.  Here  the  reduction  in  gust  sensitivity 
achieved  by  the  controller  does  not  fulfill  the  expectations. 


Fig.  16:  ID  lateral  gust  response  - lateral  acceleration  in  the 
cockpit;  open  and  closed  loop  sensitivity,  Ma=0.7, 
A=25000ft,  p-controller;  the  gust  sensitivity  in  the 
rigid-body  mode  frequency  range  is  reduced; 


Fig.  14  and  15  compare  the  open  and  the  closed-loop  cockpit 
vertical  acceleration  gust  responses.  As  expected,  a sensitivity 
reduction  can  be  observed  in  the  frequency  range  of  the  rigid 
body  motion.  Further,  the  integrated  flight  and  aeroelastic  con- 
troller achieves  a significant  reduction  of  the  acceleration  asso- 
ciated with  engine  and  fuselage  modes  over  the  entire  aircraft 
structure.  The  vertical  gust-induced  acceleration  associated 
with  e.g.  the  1st  wing  bending  mode  is  significantly  reduced  on 
the  wing  of  the  aircraft  while  no  improvement  (or  even  a dete- 
rioration) can  be  detected  in  the  cockpit. 

The  influence  of  the  controller  on  gust  sensitivity  is  essentially 
similar  for  both  the  uniform  gust  model  used  for  design  and  the 
ID  standard  gust  model,  although  the  quantitative  results  dif- 
fer. 


While  the  gust  sensitivity  in  the  frequency  range  associated 
with  the  dutch  roll  mode  is  somewhat  reduced  and  shifted  to 
higher  frequencies  (as  dutch  roll  frequency  is  increased),  the 
sensitivity  in  the  aeroelastic  range  is  reduced  only  for  some 
modes. 

Fig.  17  shows  the  detailed  results  for  both  the  uniform  and  the 
ID  gust  models.  It  can  be  noted,  that  the  significant  reduction 
of  the  sensitivity  peak  associated  with  the  fuselage  bending 
mode  in  the  case  of  uniform  gusts  is  not  achieved  for  ID  gusts. 
As  the  difference  in  phase  between  uniform  and  ID  lateral 
gusts  unfolds  over  the  length  of  the  fuselage,  it  is  not  surprising 
that  the  fuselage  bending  mode  is  the  most  affected.  Regarding 
the  engine  modes,  the  larger  sensitivity  reduction  can  be 
observed  for  the  ID  gust  model.  For  the  fish-tailing  mode,  the 
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controller  produces  no  improvement  for  uniform  and  even  a 
deterioration  for  ID  gusts. 

Therefore  further  research  should  concentrate  on  improving 
the  performance  with  respect  to  fish-tailing.  Further,  it  is 
expected  that  using  the  ID  gust  model  for  control  design  will 
improve  the  achievable  sensitivity  reduction,  especially  with 
respect  to  fuselage  bending. 


Frequency 


Fig.  17:  Open  and  closed  loop  sensitivity  to  lateral  gusts; 
a)  uniform  gust  model;  b)  ID  gust  model; 
p-controller  lateral  motion,  Ma=0.7,  /r=25000ft;  sensi- 
tivity reduction  depends  strongly  on  the  gust  model 

Gust  sensitivity  analysis  using  the  2D  gust  model 

In  the  analysis  run,  the  2D  gust  model  described  in  section  5 is 
employed  to  compare  the  open  and  closed  loop  gust  sensitivity. 
Only  the  longitudinal  controller  is  implemented. 


Fig.  18:  Open  and  closed  loop  gust  sensitivity  to  vertical  gusts; 
2D  gust  model;  longitudinal  p-controller; 

a)  PSD  for  outer  wing  vertical  acceleration; 

b)  PSD  for  outer  engine  lateral  acceleration; 

Fig.  18a  shows  the  open  and  closed  loop  outer  wing  vertical 
acceleration  PSDs.  The  outer  engine  lateral  acceleration  PSDs 
are  given  in  Fig.  18b.  The  reduction  of  the  sensitivity  peaks 
associated  with  the  symmetric  modeshapes  corresponds  to  the 
reduction  observed  for  the  ID  gust  model.  The  increased 


damping  of  the  wing  bending  and  engine  modes  entails  a gust 
sensitivity  reduction. 

Naturally,  the  longitudinal  controller  cannot  reduce  the  excita- 
tion of  antimetric  modeshapes.  Clearly,  for  the  example  air- 
craft, lateral  aeroelastic  control  action  is  required  for  vertical 
gust  sensitivity  reduction.  It  should  be  noted  however,  that  by 
considering  (uniform)  roll  gusts  for  lateral  control  design,  the 
influence  of  vertical  gusts  on  the  lateral  motion  has  already 
been  partly  accounted  for. 

10  Conclusion 

In  this  paper  the  response  of  a large  aircraft  to  gust  input  has 
been  evaluated.  For  the  analysis,  a high  fidelity  structural  and 
aerodynamic  model  that  also  represents  aerodynamic  coupling 
between  rigid  body  and  structural  modes  is  employed.  Using 
gust  models  of  different  accuracy  and  complexity  (uniform, 
ID,  2D,  vertical  and  lateral),  it  is  shown  that  aircraft  reaction 
can  diverge  significantly  for  large  aircraft. 

Aeroelastic  vibrations  of  the  aircraft  are  considerably  excited 
by  vertical  and  lateral  gusts.  Therefore,  realistic  gust  models 
that  describe  the  distribution  of  the  gust  over  the  aircraft  struc- 
ture are  required.  From  the  comparison  of  uniform  and  1 D gust 
responses  in  the  aeroelastic  frequency  range  (e.g.  fuselage 
bending),  it  can  be  concluded  that  uniform  gust  approximations 
are  inadequate  for  sensitivity  analysis.  This  is  especially  true 
for  large  airplanes,  where  the  delays  between  the  gust  effect  at 
the  front  and  the  rear  of  the  aircraft  are  more  important. 

Investigations  often  neglect  the  response  of  aircraft  due  to  lat- 
eral gusts,  which  also  lead  to  a considerable  excitation  of  the 
airplane.  In  particular,  the  bending  modes  of  the  fuselage  and 
the  wing  versus  fuselage  “fish  tailing”  mode  have  to  be  stud- 
ied, as  they  cause  high  loads  at  the  rear  fuselage. 

The  implementation  of  2D  gust  models  adds  new  aspects  to  the 
problem,  especially  the  response  of  the  wing  is  further  modi- 
fied. Vertical  gusts  excite  antimetric  modes  and  the  energy  of 
the  gust  is  distributed  on  more  modes.  This  explains  the  PSD- 
reduction  associated  with  the  first  wing  bending  mode.  On  the 
other  hand  there  are  modes  whose  response  is  amplified  by  a 
2D  gust.  In  particular,  if  the  gust  velocity  distribution  has  the 
same  form  as  certain  mode  shapes  (e.g.  higher  order  wing 
bending  or  engine  modes),  an  amplification  of  the  correspond- 
ing PSD  results. 

The  closed  loop  system  analysis  (section  8)  demonstrates  that 
the  integrated  controller  meets  the  formulated  performance 
requirements.  The  comparison  of  open  loop  and  closed  loop 
gust  induced  acceleration  PSD  curves  shows  that  aeroelastic 
control  achieves  a significant  reduction.  The  use  of  uniform 
gust  approximations  for  control  design  is  not  satisfactory,  as 
the  sensitivity  reduction  results  are  different  for  the  uniform 
and  the  ID  gust  models.  In  particular,  ID  gust  sensitivity  peaks 
associated  with  fuselage  bending  modes  cannot  be  reduced  by 
the  controller.  This  is  not  unexpected  as  the  phase  difference 
between  uniform  and  ID  gust  excitations  is  strongest  over  the 
length  of  the  fuselage. 

Therefore,  future  research  should  concentrate  on  developing  a 
realistic  time-domain  gust  formulation  that  can  be  employed 
for  state-space-model-based  control  design.  It  is  expected  that 
this  step  will  entail  a significant  improvement  in  the  controller 
performance,  especially  with  respect  to  the  fuselage  bending 
modes  and  the  fish-tailing  mode. 
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